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ABSTRACT 

We use Herschel Space Observatory data to place observational constraints on the peak and Rayleigh- Jeans slope of dust emission 
observed at 70-500 fim in the nearby spiral galaxy M81. We find that the ratios of wave bands between 160 and 500 ^m are 
primarily dependent on radius but that the ratio of 70 to 160 ^m emission shows no clear dependence on surface brightness or 
radius. These results along with analyses of the spectral energy distributions imply that the 160-500 ^m emission traces 15-30 K 
dust heated by evolved stars in the bulge and disc whereas the 70 /^m emission includes dust heated by the active galactic nucleus 
and young stars in star forming regions. 
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1. Introduction 



The Herschel Space ObservatorjQ (jPilbratt et all |2010[ ) 
provides an unprecedented view of the far-infrared and sub- 
millimetre emission from nearby galaxies. At wavelengths 
of 7 0-160 /im, the PACS instrument (iPoglitsch et al.l 
I201Q) can produce images with resolutions of 6"- 12" that 
are superior to what can be achieved with the Spitzer 
Space T elescope. At 250-500 /im, the SPIRE instrument 
(jGrifRn et al. 2010) produces images with unprecedented 
sensitivities to diffuse and point-like submillimetre emis- 
sion. We can use these data to construct spectral energy dis- 
tributions (SEDs) that sample the peak and Rayleigh- Jeans 
side of thermal dust emission, thus allowing us to probe 
the coldest dust components in nearby galaxies and place 
superior constraints on dust temperatures and masses. As 
part of the Very Nearby Galaxies Survey (VNGS), we have 
imaged the spiral galaxy M81 (NGC 3031) at 70, 160, 
250, 350, and 500 fim with PAC S and SPIRE. M81 is a 
nearby (3.63 ± 0.13 Mpc: iFreedi^ et al.l l2001[) SA(s)ab 
(Ide Vaucouleurs et al.l 19911 ) galaxy at an inclination of 
59.0° (|de Blok et al.l l2008r with well defined spiral arms. 
The Herschel data allow us to extract SEDs for ^ 0.7 kpc 
subregions that are small enough that we can distinguish 
arm and interarm regions within M81. We use these data 
to explore the dust temperatures and masses and to under- 
stand the heating sources for the dust. 



^ Herschel is an ESA space observatory with science instru- 
ments provided by Principal Investigator consortia. It is open 
for proposals for observing time from the worldwide astronomi- 
cal community. 



2. Observations and data reduction 

The PACS observations were performed as four pairs of 
orthogonal scans covering 40' x 40' using a 20" s^^ scan 
rate. PACS can perform simultaneous observations in only 
two wave bands; we chose the 70 and 160 /im bands 
since they were expected to bracket the peak of the SED 
better. The data were reduced using a combination of 
an adapted Herschel Interactive Processing Environment 
(HIPE) 3.0 pipeline and Scanamorphos (Roussel et al. in 
prep.). Starting from the raw detector timelines, HIPE was 
used to mask dead and saturated pixels, convert the signal 
to Jy pixel^^, and remove cosmic rays. Scanamorphos was 
then used to map the data while simultaneously removing 
1// drifts in the signals. Finally, we subtracted the median 
backgrounds from the images. The photometric calibration 
has an accuracy of 10% at 70 /im and 20% at 160 /im, and 
the full-width half-maxima (FWHM) of the 70 and 160 /im 
point spread fun ctions (PSFs) are 6" and 12", respectively 
(|Poglitsch et al.ll2010[ ). The RMS noise levels are 0.12 mJy 
arcsec"^ in both the 70 and 160 /im bands. 

The SPIRE observations were performed as two orthog- 
onal scans covering 40' x 40' using a 30" s~^ scan rate. A 
modified HIPE 3.0 detector timeline pipeline was used to 
remove cosmic rays, flux calibrate the data, a nd apply tem- 
perat ure drift and response corrections (see iPohlen et al.l 
I2OIOI for details). We then removed offsets between the 
detector timelines in two steps. First, we subtracted the 
median signal from each bolometer observed during the en- 
tire observation. Then we applied an iterative process to 
remove residual baseline signals that appear as stripes in 
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Fig. 1. 70-500 /im images of M81 covering 20' x 30' with north up and east to the left. The images are scaled logarith- 
mically. The green circles in the lower left corner o f each image show the FWH M of the PSF. The cyan ellipse in the 
250 ^m image shows the D25 isophote (26'. 9 x 14'.l: lde Vaucouleurs et al.]ll99l[ ): the radius is equivalent to 14 kpc. The 
blue squares in the 250 /im image show the 42" regions for which SEDs are plotted in Fig. U) 
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Fig. 2. Images of the 70/160, 160/250, 250/350, and 350/500 /im surface brightness ratios in the optical disc of M81. 
Each image is created using data with PSFs that match the 500 /J,m PSF. Regions not detected at the 3cr level in the 
two bands used for each ratio are left blank. The image sizes and orientations are the same as for Fig. [TJ The circles in 
the lower left corner of each image show the FWHM of the 500 //m PSF. 



the maps. In this process, we first created a map. Then, 
for each bolometer timeline in each scan leg, we measured 
the signal in the map that corresponded to the bolome- 
ter's position, we calculated the median difference between 
the bolometer signal and the corresponding map signal, 
and we subtracted this function from the bolometer sig- 
nal. These steps were repeated 40 times to completely re- 
move stripes from the data. Finally, we subtracted median 
background signals from the images. The resulting images 
have flux calibration uncertainties of 15%, and the 250, 350, 
and 500 ^ m PSFs have FWHM of 18", 25", and 37", re- 
spectively (jSwinvard et al.ll2010l ). The RMS noise levels are 
0.040, 0.019, and 0.008 mJy arcsec'^ in the 250, 350, and 
500 /Ltm bands, respectively. 

To create ratios of surface brightnesses measured in two 
wave bands, we matched the PSFs, which we treated as 
Gaussian, to the PSF of the 500 /im data. For statistical 
analyses on surface brightness ratios and for creating SEDs 
of subregions within the galaxies, we then rebinned the data 
in all images into 42" (~ 0.7 kpc) square pixels (selected 
because it is an integer multiple of the 500 /im pixel size 



that is larger than the PSF FWHM for the 500 ^m data). 
For these analyses, we only used used 42" pixels with ia 
detections in all bands. 



3. Results 

Figure [1] shows the structures traced by the various 
Herschel wave bands, which look very similar to each other 
and to the 5.7-24 ^m Spitzer images (jGordon et al.ll200l : 
IWillner et al.l[2004l) . All images trace the same spiral struc- 
ture and individual infrared sources within the disc of the 
galaxy. Diffuse, extended sources detected outside the op- 
tical disc of M81 in the Herschel data seem most likely to 
be associated with dust in the Milky Way (Davies et al. in 
prep.). 

To understand the heating mechanism for the dust, we 
examined how surface brightness ratios varied with surface 
brightness and with radius. Variations with surface bright- 
ness would suggest that the dust is heated locally and that 
the emission is linked to star formation, whereas radial vari- 
ations in the ratios would indicate that the dust emission 
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Fig. 3. The 70/160, 160/250, 250/350, and 350/500 fim 
surface brightness ratios versus the surface brightness (left 
and center) and incUnation-corrected galactocentric radius 
(right). The data were measured in 42" subregions in im- 
ages with PSFs that matched the PSF of the 500 /im data. 
Best fit hues are shown for all plotted data except for two 
relations involving the 70/160 fiui ratio, where the fits were 
very poor; corresponding slopes are given in the panels. 
The R values are the Pearson correlation coefficients for 
the plotted data. Note that, in the left-side and center pan- 
els, the logarithm of the surface brightnesses are used for 
the best fit lines and correlation coefficients. 



is more strongly affected by the evolved stellar populations 
in the bulge and disc. Figure [2] shows images of the 70/160, 
160/250, 250/350, and 350/500 /im surface brightness ra- 
tios. Additionally, Fig. [3] shows how the ratios measured 
in 42" (~ 0.7 kpc) subregions vary with surface brightness 
and galactocentric radius. 

The absolute value of the correlation coefficient R 
for the relations between radius and either the 160/250, 
250/350, or 350/500 /im ratios is generally higher than that 
for the relations between surface brightness and the ratios, 
which shows that these ratios are more strongly depen- 
dent on radius (although the 160/250 /im ratio may also 
be partly dependent on 160 /im surface brightness based 
on the high value of R). This is consistent with the weak 
or absent infrared-bright regions or spiral structure in the 



1000 






Global ; 


100 


^ ■/ 


/^ 


X T=17.8 
N ±0.6 K 




^ y 




\ : 


in 


/ 






10.U 
1.0 


: ,■ 


"\ 


Nucleus ! 

. T=28.9 - 
\ ±0.8 K 


?01 








53.0 




\ 


Arm Region 


1.0 


r / 




\ ifo°r 


0.3 


;,/„ 




\ ; 


1.0 


: / 


Interarm Region - 

/ t \, T=17.4 '- 
/ \ ±1.4 K: 




/ 




^ 


0.1 


-,,,', 




, ■♦,,.,,. 











1.5 


h + 




Global - 








1.0 


-_ + 




^ 


0.5 


-- 




-_ 


0.0 


^ + 


1 


M + - 


















_ 0.4 


- 




Nucleus _ 


>^ 






- 


^ 0.2 


_ 




- 








1 


giO.O 


- 1 


' 




- 


i-0.2 








_ 




r"l" 






^ 0.5 




Arm Region - 








Eu 






- 


g 






1 ; 


f 0.0 


- 


t 


1 


t : 


1.0 


- + 


Interarm Region _ 


0.5 


- 




- 


0.0 


- 


t 


U,L„r 



100 



1000 



100 



1000 



\ (|im) 



?.(Hm) 



Fig. 4. On the left, the global SED as well as the SEDs 
for the three 42" {^ 0.7 kpc) regions shown in Fig.[T] The 
SEDs for the subregions were measured in data with PSFs 
that matched to the PSF of the 500 /^m data. The grey line 
is the blackbody modified with a A~^ emissivity function 
fit to the data. On the right are the residuals from the fit 
in logarithm space. 



images of the 160/250, 250/350, and 350/500 /im ratios. 
Moreover, the R? values, which equal the fraction of the 
variance in the data that can be accounted for by the best 
fit line, indicate that >70% of the variance in the 160/250, 
250/350, and 350/500 /im ratios can be accounted for by 
the relation with radius. 

In contrast, the 70/160 /im ratio does not vary monoton- 
ically with radius except within 2 kpc, a region in which the 
gradient in the 160/250 /im ratio also increases. This could 
represent enhanced dust heating within this radius that is 
powered by the active galactic nucleus (AGN), by strong 
central star formation activity, or by the bulge stars, which 
have a high central density. The 70/160 /im ratio versus 
160 /im surface brightness exhibits no obvious trend and 
only a statistically weak trend (with i?^ < 0.3) is visible 
in the plot of the 70/160 /im ratio versus 70 /im surface 
brightness, although the best fit line poorly describes the 
data. None the less. Fig. [5] shows that the ratio increases to 
>0.3 in infrared-bright regions in the spiral arms. 

Figure |3] shows the SED integrated across the optical 
disc (with supplement al 60 and 100 /im IRAS data added 
from JRice erall (|1988[ )) as well as the SEDs for the 42" re- 
gions centered on the nucleus, and examples of an infrared- 
bright source and an interarm region. Based on visually in- 
specting and fitting functions to the SEDs, these example 
regions were typical to similar regions at similar radii. In the 
nucleus, we were able to fit a single blac kbody modified with 
a A~^ emissivity function (based on the lLi fc Drainej (|200lD 
emissivity function) to the 70-350 /im data, but the 500 /im 
data point could not be fit with the same thermal compo- 
nent, although the mismatch between the fit and model is 
only 2(7. This result and the low 350/500 /im ratio for the 
nucleus seen in Fig. |3] (which is 3cr below the best fit line) 
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suggest that the 500 jjlth nuclear emission hkely includes a 
non-thermal component associated with the AGN. Based 
on the SED fit, we estimate that the non-thermal 500 /xm 
emission is 0.05 ± 0.03 Jy, which is ~ 2.5 x below a power 
law extrapola t ed fro m the mm and cm data presented by 
iMarkoff et al.l ()2008[ ). The discrepancy could be explained 
by the low signal-to-noise in the estimate from the SED 
fit or by variability in the AGN emission; th e 870 /im flux 
densi ty has been observed to vary by 3x ([Markoff et al.l 
1200811 . 

In the other SEDs, we found that single blackbodies 
modified with X~^ emissivity functions could be fit accu- 
rately to the >100 /im data without the fit overpredicting 
the observed 70 /im measurement, but fits that included 
the 70 /^m data point did not accurately replicate the peak 
of the SED. No evidence is found for the excess emis- 
sion at submillimetre wavelengths sometimes attributed to 
dust with < 10 K temperatures or shallow emissivities, al- 
though prior results had indicated tha t this emission would 
be more prominent at >500 /im (e.g. Galliano et al. '2005!; 
LBendo et al. 2006; Zhu et al. 2009; O'Halloran et al. 2010). 
By applying to the data for the global SED the equation 
Mdust = \fvD'^M\K B{v,T)\ (wh e re D is distance, Ki, is the 
dust opacity from iLi fc Draind ()2001[ ). and B{v,T) is the 
best fitting modified blackbody), we estimated the global 
dust mass to be 3.4 ± 0.5 x 10^ M(7i . Given that the at omic 
gas mass is 3.64 ± 0.18 x 10^ Mq (|Walter et al.ll2008l) and 
the m olecular gas mass is negligible in comparison (jSagd 
Il993l Sanchez-Gallego et al. in prep.), we estimate that the 
gas-to-dust ratio is 107 ± 17, which is within the range of 
'^ 100-200 expected for solar metallicity objects based on 
the depletion of metals from the gaseous phase of the in- 
terstellar medium or comparisons of gas column density to 
dust extinction (e.g Whittct 2003). Hence, this simplistic 
modified blackbody fit may be a fair representation of the 
emission from the bulk of the dust mass in M81, although 
more sophisticated modeling should not only yield more ac- 
curate masses but also describe the emission from warmer 
dust components. 

The SED fits along with the results from Figs. [Hand [3] 
imply that the 70 /xm band traces dust heated by a differ- 
ent source than the dust that primarily emits in the 160- 
500 /im bands. Although the 70/160 /im ratio exhibits a lot 
of scatter, the enhancements in the 70/160 /im ratio in the 
spiral arms implies that the 70 /im band may be affected by 
star formation on local scales. Meanwhile, the radial vari- 
ations in the 160-500 /im ratios and the SED fits suggest 
that ~ 20% of the 60 /im emission, ~ 30% of the 70 /im 
emission, and ~ 100% of the > 100 /im emission originates 
from dust heated by evolved disc and bulge stars. This is 
consistent with prior results suggesting that ~ 5-100% of 
the 60 and 100 /zm emission from nearby galaxies originates 
from dust heated by evolved stars ( e.g. ISauvaee &: ThuanI 
ll992tlWalterbos fc Greenawaltll 19961 ). If this interpretation 
is correct, we anticipate that dust emitting at 160-500 /im 
in other galaxies with relatively large fractions of old stars 
(E-Sab galaxies) will also have 160-500 /xm colours that 
depend upon radius, but galaxies with relatively large frac- 
tions of young stars (Sc-Im galaxies) will have 160-500 /xm 
colours that may depend more on infrared surface bright- 
ness, as heating by the evolved stellar population becomes 
insignificant. The results also imply that the conversion of 
infrared fluxes integrated over very broad rang es (e.g. 8- 
1000 /xm) to star formation rates, as done by IZhu et all 



(|2008D . lRieke et al.l([2009h . and lKennicutt et al.l(|2009[ ). wih 
be accurate as long as the integrals contain a significant 
amount of emission shortward of 160 /im that traces dust 
heated by star formation. However, it may not be possible 
to derive accurate star formation rates from dust emission 
measured solely at >160 /xm. 

In conclusion, these results for M81 demonstrate how 
Herschel 70-500 /im data can be used to not only measure 
more accurate dust temperatures and masses but also de- 
termine the dust heating sources in nearby galaxies. Further 
work with data from the VNGS and other surveys will al- 
low us to determine whether dust traced by the 160-500 /im 
bands in other spiral galaxies is also heated by evolved 
stellar populations and whether variations in the relative 
strength of dust heating by evolved stars varies across the 
Hubble sequence. 
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